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Abstract 

We report on the first systematic study of spin transport in bilayer graphene (BLG) as a function 
of mobility, minimum conductivity, charge density and temperature. The spin relaxation time t s 
scales inversely with the mobility \i of BLG samples both at room temperature (RT) and at low 
temperature (LT). This indicates the importance of D'yakonov - Perel' spin scattering in BLG. Spin 
relaxation times of up to 2 ns at RT are observed in samples with the lowest mobility. These times 
are an order of magnitude longer than any values previously reported for single layer graphene 
(SLG). We discuss the role of intrinsic and extrinsic factors that could lead to the dominance of 
D'yakonov-Perel' spin scattering in BLG. In comparison to SLG, significant changes in the carrier 
density dependence of r s are observed as a function of temperature. 
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FIG. 1: (Color online) (a) AFM image of a BLG sample after MgO deposition: rms roughness ~ 
0.3 nm. (b) SEM image of a BLG sample with multiple non-local spin valves. 

The demonstration of micrometer long spin relaxation lengths in graphene by Tombros 
et al. pQ has made this two-dimensional material an extremely promising candidate for 
spintronics applications. So far most spin transport studies have focused on single layer 
graphene (SLG) [THP3] while the equally important bilayer graphene (BLG) has not yet 
received much attention. This is surprising since BLG has unique electronic properties 
which differ greatly from those of SLG (effective mass of carriers, electric-field induced band 
gap) and also differ from those of regular 2D electron gases (chirality) [TH IT5] . It is currently 
believed that spin relaxation in SLG is limited by the momentum scattering from extrinsic 
impurities [21 [5j [16] . Unlike SLG, the scattering from such charged impurities is reduced by 
the enhanced screening in BLG [T7]. This leads to a relative importance of short-range (SR) 
scatterers in determining the transport properties such as the temperature ( T) and charge 
carrier density (n) dependence of the BLG conductivity (cr) [TTHT9] . In addition, interlayer 
hopping also plays an important role in the electronic properties of BLG and is predicted 
to cause an enhanced intrinsic spin-orbit (SO) coupling (up to 0.1 meV in clean samples) 
in comparison to SLG [20J. Since charge and spin transport are highly linked, it is natural 
to expect a profound difference in the nature of spin transport in BLG. In particular, the 
unique electronic properties of BLG may offer new avenues to manipulate the spin degree 
of freedom. 

In this Letter we report on spin transport studies in BLG both at room temperature 
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FIG. 2: (Color online) RT data: (a) a vs. Vq and t s vs. Vq for BLG. (b) Non-local resistance 
as a function of the in-plane magnetic field B y (T). The blue and red arrows show the field sweep 
direction while the black arrows show the relative magnetization orientations of the injector and 
detector electrodes. Hanle precession measurement for a perpendicular magnetic field B Z (T) sweep 
for (c) the same sample with fi ~ 2000cm 2 /Vs and (d) for a sample with fi ~ 300cm 2 /Vs. 

(RT) and at low temperature (LT) using MgO barriers. Spin valve devices in the non-local 
geometry are fabricated on two types of MgO-covered exfoliated graphene samples using 
standard e-beam lithography techniques. For global MgO samples, MgO covers the entire 
graphene surface while for local MgO samples the MgO is only under the Co electrodes. 
This is followed by the evaporation of the ferromagnetic contacts (Co) [2T]. The AFM 
image after MgO deposition and the SEM image after the device fabrication for one of the 
BLG samples are shown in Fig.l. In order to investigate the nature of spin scattering in 
BLG, we have evaluated the spin relaxation time t s as a function of four paramters: (1) the 
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field-effect mobility fi, (2) the minimum conductivity a m in, (3) the charge carrier density 
n and (4) the temperature T in the range of 5 K to 300 K. Among these parameters the 
mobility dependence of t s provides the most direct way to deduce the dominant scattering 
mechanism: a linear dependence of t s on \i (or r p ) is a priori suggestive of an Elliott- Yafet 
(EY) spin scattering mechanism [22J, while the inverse relation (r s oc l/fi oc l/r p ) will 
indicate the dominance of D'yakonov-Perel' (DP) like spin scattering mechanisms [23J. In 
general, both mechanisms could be simultaneously relevant. For this study, we have selected 
representative 17 devices on 6 BLG samples whose field effect mobilities vary by more than 
one order of magnitude, from fi ~ 200 cm 2 /Vs to 8000 cm 2 /Vs. The effect of impurity 
scattering on r s can also be deduced from a m i n , which depends on impurity concentration. 
The third parameter chosen is the charge carrier density: the density dependence of r s and 
t p is used to identify the spin scattering mechanism in SLG (EY). In BLG, the r p is often 
taken to be a constant under the assumption of charge scattering from weak short-range 
scatterers and charged impurities. Recent LT experiments j2H EH] indicate, however, that 
the density dependent r p results from strong short-range scattering in BLG. As a general 
trend in most of our samples a is linear at high carrier density 1-3 x 10 12 /cm 2 (see Fig. 2 
upper panel). As the variation of r p is weak in this charge density range, we particularly 
compare the product t s t p and the ratio t s /t p as a function of carrier density n; the former 
is expected to be a constant for DP whereas the latter is a constant for EY. The fourth 
parameter chosen is temperature: unlike SLG, there is a strong temperature dependence of 
charge transport in BLG, which should also reflect itself in spin transport parameters. 

Prior to any spin transport measurements, we have characterized the BLG device conduc- 
tivity as a function of back gate voltage (Fig. 2a) using the standard four terminal geometry. 
As also observed by other groups [H [10], the graphene samples are electron doped in the 
spin- valve configuration due to doping by the Co/MgO barrier [26]. The non-local spin 
signal (AR) is measured by sweeping the in-plane magnetic field in a loop from negative 
(-80 mT) to positive (80 mT) and then back to negative values (-80 mT) [27]. A clear 
bipolar spin transport signal is observed at RT (see Fig. 2b), with a positive value of the 
non-local resistance for parallel alignment of the electrodes' magnetization and a negative 
resistance for the anti-parallel alignment (AR = 4 Q). To confirm the observed spin signals, 
conventional Hanle spin precession measurements [2H] are performed at the same electron 
density of n = 1.5 x 10 12 /cm 2 above which the conductance is linear and the mobility is 
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FIG. 3: (Color online) Results of Hanle precession measurements at RT for 17 BLG spin devices 
with mobility varying from 200 - 8000 cm 2 /Vs. The data points with the same symbol represents 
different junctions on the same sample. Open and closed symbols correspond to samples with 
global and local MgO, respectively, (a) t s taken at n = 1.5 x 10 12 /cm 2 vs. \i plotted on a log-log 
scale, (b) t s taken at the CNP vs. a m i n for BLG samples at RT. 

well defined within the Boltzmann approximation (Figs. 2c and 2d). The magnetic field 
dependence of the non-local resistance is fitted by 

f°° 1 —iP 1 -t 

R nI (x / — ===e«^ cos(cu L t)e^dt (1) 

where D s is the spin diffusion coefficient and ul = gfiB^/h is the Larmor frequency for spin 
in an external magnetic field and thus gives the values for the spin parameters |28j . At n = 
1.5 x 10 12 /cm 2 we obtain a diffusion constant of D s = 0.0032 m 2 /s and a spin relaxation 
time of r s = 135 ps for a sample with /x = Acr/eAn = 2000 cm 2 /Vs (Fig. 2c). The values 
give a spin relaxation length of 0.7 /jm. The spin relaxation time t s as a function of gate 
voltage (doping) is plotted in Fig. 2a (lower panel) showing an increase (< 10%) of t s with 
doping, away from the charge neutrality point (CNP). This is qualitatively similar to the 
gate tunability of t s in SLG, although r s shows a weaker dependence in BLG at RT. 

We next evaluate in detail the dependence of t s on \i in the Boltzmann regime. Such 
mobility dependent studies have been useful in identifying spin scattering mechanisms in 
inorganic semiconductor systems |29j. As shown in Fig. 3a by a log-log plot, we observe 
an inverse dependence of r s on the mobility. Note that all data are taken at an electron 
density n = 1.5 x 10 12 /cm 2 . In samples with the highest mobility t s is only 30 ps. On 
the other hand, we observe a spin relaxation time t s of up to 2 ns at RT for samples 
with the lowest mobilities (see corresponding Hanle curve in Fig. 2d, [30J). Such values 
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for t s are one order of magnitude longer compared to values reported so far in any SLG 
experiment. Furthermore, this strong variation of r s with fi offers the most direct evidence 
of the correlation between spin and charge transport. Since higher mobilty samples will 
typically involve higher momentum relaxation time t p , assuming \x oc r v in the Boltzmann 
regime, the inverse dependence of r s on \i clearly demonstrates that the DP mechanism is 
the dominant spin scattering mechanism in BLG at RT. 

The same plot of r s vs. fi is also useful to elucidate the possible origin of the DP 
mechanism in our samples. For this, we first estimate the strength of spin-orbit (SO) coupling 
(A) in BLG using the expression: l/r s = ^e// r p = 4A 2 r p //i 2 , where fi e // is the effective 
Larmor frequency of the precessing spins and A is the corresponding SO coupling strength 
[16] . The value fi e // is obtained from the r s vs. /i data (Fig. 3a), Q e ff — 407 ± 25 GHz, 
which gives A ~ 0.14 ± 0.01 meV. Additionally, A only weakly depends on temperature (21] . 
Therefore, it is unlikely that low energy phonons (such as acoustic phonons) are responsible 
for the observed spin scattering. The most important open question is whether this SO 
coupling is intrinsic or extrinsic in nature. The intrinsic SO coupling of BLG is expected to 
lead to Aintrinsic up to 0.1 meV in clean BLG samples [20]. While this is in good agreement 
with the value extracted from the t s vs. /i plot, the influence of both the externally applied 
electric field and the role of adatoms cannot be excluded. Since interlayer hopping is involved 
in BLG, electric field dependent modifications to the intrinsic SO coupling are expected. 
Adatoms, on the other hand, induce local curvature to an otherwise flat graphene lattice 
and can cause spin scattering by both EY and DP mechanism [T21 [TB]. In SLG, the recent 
studies on the influence of external adatoms on spin transport also shows an increase in t s 
with increase in adatom concentration indicating a DP like scattering at LT [9]. However, 
in the case of BLG the role of the adatoms, in determining DP or EY spin scattering, might 
be even smaller due to a higher lattice stiffness [31], thus reducing the adatom induced 
SO coupling strength. Moreover, the effect of the charged impurities is also reduced due 
to enhanced screening in BLG. Thus the two prominent factors (charged impurities and 
adatoms) responsible for the spin scattering in SLG play a minor role in BLG which could 
also be the reason why we see longer spin relaxation times in BLG. More investigations are 
needed to clarify the type and concentration of spin scatterers in single and bilayer graphene 
and to differentiate the contribution from extrinsic and intrinsic factors to the SO coupling 
in BLG samples with disorder. 
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FIG. 4: (Color online) (a) r s vs. /i for 8 BLG junctions at 5K. (b) Upper panel: R vs. n for BLG; 
Lower panel: r s vs. n for T = 300K (black circles), 50K (blue circles) and 5K (red circles), (c) 
Upper panel: R vs. n for SLG; Lower panel: t s vs. n for SLG at RT (black circles) and at 5K (red 
circles), (d) Upper panel: t s vs. T for four densities , n = CNP (black circles), 0.7 (red circles), 
1.5 (blue circles) and 2.2 (green circles) x 10 12 /cm 2 . Lower panel: t s vs. Tfor SLG at n = CNP 
(black circles), 1.5 (blue circles) x 10 12 /cm 2 . 

We next evaluate the dependence of r s on the RT minimum conductivity a m i n and charge 
carrier density n of the BLG samples. The r s vs. o m in graph (Fig. 3b) indicates a decrease 
in t s with increasing a m i n . The value of u m in a t RT varies significantly between clean and 
dirty samples and is higher for cleaner samples [121 E21 E3] ■ The inverse dependence of t s on 
u m in indicates a higher spin relaxation time in dirtier samples. Therefore, this correlation is 
in good agreement with the above conclusion that DP spin scattering is dominant in BLG at 
RT. Finally, we analyze the weak dependence of r s on n for individual samples (Fig. 2a). At 
first we note that at RT the density dependence of t p = am*/ne 2 shows a gradual decrease 
with increasing n in the density range 1 - 3 x 10 12 /cm 2 for our samples [21]. In this range, 
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the quantity t s t p is almost constant (~ 4% change with charge density) while t s /t p shows an 
increase of 28% with increasing charge density at RT. This is consistent with the dominance 
of the DP mechanism at RT. 

So far we have evaluated spin transport at RT. Next we evaluate spin transport as a 
function of temperature. Here, it is important to note that for BLG (Fig. 4b), and unlike SLG 
(Fig. 4c), the charge transport already shows strong changes with decreasing temperature 
due to the thermally activated nature of carriers near the CNP [211 Ell [35] . We may expect 
these changes to be reflected in the spin transport, as the temperature is lowered. We 
first note that for charge transport, there are two distinct bilayer specific density regimes 
(Fig. 4(b) upper panel): (A) n > 1 x 10 12 /cm 2 , involving temperature independent transport 
and (B) n < 1 x 10 12 /cm 2 , involving thermally-activated transport [21]. For regime (A) 
away from the CNP the mobility is well defined. Here the t s can be easily evaluated, similar 
to the approach used for analyzing RT data. We measured 8 devices in 3 samples with /i 
ranging from 700 cm 2 /Vs to 3800 cm 2 /Vs at T = 5 K. The data plotted in Fig. 4a, show 
that the inverse dependence of t s on /i persists down to 5 K, demonstrating DP as the 
dominant spin scattering mechanism even at LT. At the same time, we note that the density 
dependence of r s and t p is rather weak in the Boltzmann regime for n > 1 x 10 12 /cm 2 
and does not allow for a clear assignment of the dominant spin scattering mechanism. The 
quantities t s t p and r s /r p show comparable small changes in this regime both at 50 K and 
at 5 K [21]. Therefore, in this density regime we are left only with the mobility dependence 
at LT. 

We finally consider the low density regime (n < 1 x 10 12 /cm 2 ) around the CNP, involving 
the thermally activated behavior of the charge transport. Here, the minimum conductivity 
at CNP (a m i n ) is not any more a suitable parameter for scaling t s . The available a m i n data 
scatters significantly at LT (not shown), but the overall variation in a m i n is small compared 
to RT. This is not surprising since a min is expected to take a disorder independent value 
of 3e 2 /7rh at T = K [36]. Thus we are left only with the density dependence of t s to 
elucidate on what happens near the CNP at LT. This density dependence of t s at RT shows 
a minimum at CNP, as discussed above. As the temperature is lowered, this minimum in r s 
is gradually suppressed and the slope of r s (n) changes sign. Finally, at T = 5 K the density 
dependence of t s shows strong enhancement (> 50 %) near the CNP (Fig. 4b, lower panel). 

The key to understanding which spin scattering mechanism dominates at LT near CNP 
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lies in the density dependence of r p . We note that r p = am* /ne 2 (Boltzmann approximation), 
estimated for our samples, gives a quantitative estimate only in the high density (Boltzmann) 
regime. Closer to the CNP, r p extracted from the above assumption shows an increase with 
decreasing density |21j. Recent detailed experiments on BLG at LT have shown that a 
divergence of r v near the CNP is indeed observed [211 EU EH] . With this, a correlation of 
T s [n) with T p (n) suggests a transition from DP to EY spin scattering mechanism at LT, 
around the CNP (21]. Since the momentum scattering mechanism is different near CNP at 
LT, it is not surprising that the spin scattering mechanism is also different. We note that 
two mysteries remain to be resolved: (1) why is a transition to the different spin scattering 
mechanism observed in BLG (1) near CNP, and (2) only at LT. One possible explanation 
could be related to the thermally activated nature of carriers in this density regime. 

In conclusion, we have demonstrated spin injection and detection in BLG across MgO 
barriers and observed spin relaxation times up to 2 ns at room temperature. Our systematic 
study shows that at RT spin scattering in BLG follows an inverse dependence of r s on both 
the mobilty /i and the room temperature <J m i n , indicating a DP spin scattering mechanism. 
We disuss the role of intrinsic and extrinsic factors that could lead to the dominance of 
DP spin scattering in BLG. While the inverse scaling of mobility with r s persists down to 
5 K, the density dependence of t s indicates deviations from DP mechanism at these low 
temperatures near the charge neutrality point. 

The Aachen team acknowledges support by DFG through FOR 912. B.B and M.P. 
acknowledge support from JARA-FIT. The Singapore team acknowledges support by the 
Singapore National Research Foundation under NRF RF Award No. NRF-RF2008-07, NRF- 
R-143-000-360-281, US Office of Naval Research (ONR and ONR Global), and by NUS 
NanoCore. 

Note added: After the submission of this manuscript, related works on BLG and on 
multilayer graphene became available [37, 38J. 
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I. BLG spin valve device Fabrication: 

The graphene samples for our devices are prepared by micromechanical exfoliation of 
graphite onto a Si/SiC>2 substrate. Both SLG and BLG samples are first identified by their optical 
contrast which is then confirmed by Raman spectroscopy (laser wavelength = 514 nm) (see Fig 
Sib). For global MgO samples (MgO covered over the entire graphene surface), the next step is 
to deposit MgO under molecular beam epitaxy (MBE) conditions followed by standard e-beam 
lithography for patterning the sample. Co electrodes are then deposited by MBE. For local MgO 
samples (MgO under the Co electrodes only), we use first e-beam lithography for electrode 
patterning. This is followed by a single run deposition of MgO/Co under identical MBE 
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conditions. A thin film of MgO (1-2 nm) is deposited between Co and graphene to overcome 
the conductivity mismatch and to increase the injection efficiency of spin polarized carriers. 
Prior to and after MgO deposition (deposition rate of 0.007 nm/s), the samples were annealed for 
one hour at 200°C in vacuum (5 x 10" 10 mbar). The pre- and post-annealing steps allow a uniform 
MgO growth with typical thickness of 2 nm and roughness of around ~ 0.3 nm on the graphene 
surface (see Fig. la). Although we do not use TiO [SI] as a buffer layer for uniform MgO 
growth, we do obtain a uniform MgO layer but with pin holes in some of our MgO barriers. 
Standard e-beam lithography is performed to write the electrode patterns. The widths of the FM 
electrodes in our devices are between 300 nm to 1 \im with the separation between the electrodes 
in the range 0.50 jam to 2.6 jam. In a single run we deposit a 30 - 40 nm thick Co layer for the 
electrodes and contact pads. 

It is also important to note that after MgO deposition, we do not see any amorphization of 
graphene as reported by other groups [S2]. In Ref. [S2], the authors use sputter deposition for 
A1 2 3 and MgO layers. In contrast, we have performed MBE growth of MgO. Although we 
cannot directly compare both techniques, we expect that the graphene layer structure is less 
affected during the MBE growth of MgO. This is confirmed by Raman spectroscopy. Figure Sib 
shows one such Raman spectrum with the characteristic G and 2D peaks of graphene after MgO 
deposition confirming the quality of graphene samples even after MgO deposition. 
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Figure SI: (a) Schematic of a graphene based non-local spin valve. Measurements are performed 
with standard a.c. lock-in techniques at low frequencies with currents in the range of l-10uA. (b) 
Typical Raman spectrum of a MgO covered bilayer sample with G peak at 1579 cm" 1 and 2D 
peak at 2695 cm" 1 with FWHM -55. 

II. Temperature Dependence of BLG resistivity. 

In BLG, the conductivity minimum at the charge neutrality point (CNP) shows a distinct 
behaviour when compared to SLG (see fig. S2). The reason for such an increase in the device 
resistivity with lowering temperature is attributed to the thermally activated nature of charge 
carriers near CNP [S3, S4]. This can be confirmed by analyzing the temperature dependence of 
conductivity a at different densities. Figure S2 (b) shows the temperature dependence of a at two 
different densities n (CNP, 1 x 10 12 /cm 2 ), indicating an insulating behaviour at low T [S2]. 
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Figure S2: (a) Carrier density dependence of BLG resistivity for the temperature range 5 - 300 



K. (b) Temperature dependence of conductivity at the CNP and at n = 1.5 x 10 12 /cm 



III. Temperature Dependence of t £ and T n for BLG: 

Determining the relation between the spin relaxation time and momentum relaxation time is 
essential to deduce the dominant spin scattering mechanism: a linear relation, x s a x p indicates 
the dominance of Elliott- Yafet (EY) spin scattering, while the reciprocal relation suggests the 
dominance D'yakonov-Perel' (DP) like spin scattering. Hence we plot the product x s x p and the 
ratio x s /xp as a function of density at different temperatures [Fig S3 (b, c, & d)]. Fig. S3 (a) shows 
the variation of x p as a function of carrier density. The estimate of x p assumes Boltzmann 
transport, x p = am*/ne 2 (where m* = 0.03m e for BLG). The strong increase of x p near the CNP is 
qualitatively in good agreement with the results obtained from recent experiments on BLG [S5, 
S6]. This increase has been attributed to logarithmic corrections to x p arising from the presence 
of strong short-range scatterers. Therefore, it is reasonable to assume a similar scaling of x p (n) in 
our samples at LT near the CNP. The values for x s x p at 300 K as function of carrier density in 
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Fig. S3b are almost constant (4% change with carrier density). This confirms the observed DP 
scattering. On the other hand, the ratio x s /x p shows a significant increase of 28% at room 
temperature. However, at LT both the quantities vary by a comparable amount, i.e. both r s r p and 
t s /t p show in the density range 1-3x10* /cm a decrease by ~ 10% at 50 K and by ~ 25% at 5 
K. This makes it difficult to deduce the dominant scattering mechanism from density dependence 
alone. Therefore, at LT the scattering mechanism is determined only by the mobility dependence 
of the spin relaxation time. Moreover, it should be noted that the dependencies of r s r p and xjx v 
given here are consistently observed in all our samples. 

Effect of electron-hole puddles at CNP in BLG at LT : As discussed, the logarithmic corrections 
arising from the presence of strong short-range scatterers, causes a strong increase in r p near the 
CNP in BLG. Thus, a comparison of the r p and r s values indicates a transition from DP to EY 
like spin scattering in BLG at LT near the CNP. Moreover, the presence of electron-hole puddles 
for these bias voltages may also suggest the possibility of spin scattering via the Bir-Aronov- 
Pikus (BAP) mechanism [S7]. If such electron-hole scattering together with spin flip is present, 
the spin scattering is expected to happen only at the boundaries of electron-hole puddles and not 
uniformly across the samples. This is because the BAP spin scattering requires a significant 
overlap in the electron and hole wave functions. Although we cannot completely rule out spin 
scattering of BAP -type, we believe that this is not likely because it is usually only relevant in 
hole-doped systems with large effective mass of the charge carriers [S8]. However, precise 
theoretical studies on the electron-hole scattering near the CNP are required to fully estimate the 
influence of electron-hole exchange interaction in spin scattering. 
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Figure S3: (a) The momentum relaxation time, r p = om*/ne calculated in the Boltzmann 
framework as a function of carrier density n for BLG at T = 300 K (black circles), 50 K (red 
circles) and 5 K (blue circles); (b, c and d) The carrier density dependence of the product x s x p and 
the ratio t s /t p , which identifies the dominant scattering mechanism for T = 300 K, 50 K and 5 K, 
respectively. A constant value for x s x p indicates DP while a constant value for x s /x p indicates EY 
mechanism. The arrow in fig b shows the significant change in the ratio x s /x p with density at RT 
when compared to the change in x s x p . 



IV. Determination of the strength of spin-orbit coupling as a function 
of Mobility and temperature. 
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The spin-orbit coupling (A) is estimated from the Larmor frequency as A = £2 e ffh72 [S6]. Here the 
Larmor frequency is defined as £2 e ff 2 = l/x s x p and is evaluated from fig. 3a & 3b (see main text). 
The values of A thus obtained show only weak temperature dependence (0.14 ± 0.01 meV at RT 
and 0.12 ± 0.009 meV at LT) and mobility dependence. This eliminates the possible contribution 
from low energy phonons and suggests contributions from adatoms and from intrinsic spin-orbit 
coupling to the observed spin scattering in BLG. Furthermore, the SO coupling arising from 
phonons (acoustic and surface-optical phonons) and the charged impurities in Si0 2 have been 
estimated from first principle calculations for SLG. Both these extrinsic factors were shown to 
cause SO coupling weaker by orders of magnitude [S9] than required for the observation of sub- 
nanosecond spin scattering times in SLG. 
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